Velocity and Energy Profiles In Two- vs. Three-Dimensional Channels: 
Effects of Inverse vs. Direct Energy Cascade 
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In light of some recent experiments on quasi two-dimensional (2D) turbulent channel flow we 
provide here a model of the ideal case, for the sake of comparison. The ideal 2D channel flow differs 
from its 3D counterpart by having a second quadratic conserved variable in addition to the energy, 
and the latter has an inverse rather than a direct cascade. The resulting qualitative differences in 
profiles of velocity V and energy K aa a. function of the distance from the wall are highlighted and 
explained. The most glaring difference is that the 2D channel is much more energetic, with K in wall 
units increasing logarithmically with the Reynolds number Re,- instead of being Rer-independent 
in 3D channels. 



Experimental realizations of two-dimensional (2D) tur- 
bulence are not easy to come by, but there is a tradition, 
starting with Couder and coworkers [H, to achieve 
such realizations using soap- films 0]. A number of el- 
egant experiments, starting with Goldburg and cowork- 
ers, 0, [a H) on forced soap films bounded by straight 
wires, ignited an interest in two dimensional turbulence 
in a channel geometry. Indeed, a number of simulations 
[3, and models iTl were presented to compare with 

the experimental findings. While it is understood that 
such experiments suffer from three-dimensional (3D) ef- 
fects like film thickness fluctuations [6| and friction be- 
tween the film and the surrounding air it is clearly 
worthwhile to develop a reasonable theoretical model of 
ideal 2D channel flows to be able to gauge the degree of 
closeness of experiments and theory. It is quite surprising 
that not enough had been done in determining what are 
the expected velocity and energy profiles in such ideal 
channel flows, in parallel to the very well studied 3D 
channels. The aim of this Letter is to close this gap. 

We consider stationary fully developed turbulent flow 
of a fluid of unit density in infinitely long (in the stream- 
wise direction x) 2D and 3D (infinitely wide in z direc- 
tion) channels of width 2L (in the y direction), driven by 
a pressure gradient p' — —dp/dx. The velocity field is 
denoted as U(r,i) = V{y)x + u{r,t), where V{y) is the 
mean velocity and u{r,t) the turbulent velocity fluctu- 
ations. In such geometries V{y) and all the other mean 
quantities depend only on the distance from the wall y. 
We will be interested in the profiles of one-point averages, 
with the mean shear S{y) — dV/dy, the Reynolds stress 
Wly) = —{uxUy) and the mean turbulent kinetic energy 
K{y) = (|u^|)/2 being the primary ones. In developing 
a model for these profiles it is important to maintain as 
many exact relations as possible, and one such exact re- 
lation is the momentum balance equation [l2] which is a 
direct consequence of the Navier-Stokes equations: 



vS{y) + W{y)=p'{L^y), 
in which v is the kinematic viscosity. 



(1) 



The second exact relation is the turbulent kinetic en- 
ergy balance: V ^V + e. Here V{y) = WS is the energy 
production, e{y) — v (^{djUi)^^ is the viscous dissipation 
and 'D{y) — T)^ + I?t is the energy spatial transfer, con- 
sisting of the viscous and turbulent contributions: 



X>v = -V- 
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dy^ 
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where p is the pressure fluctuation. This is as far as once 
can go exactly. Now we need to model V in terms of 
the one-point averages. This step is identical in 3D and 
2D channels; in its simplest version it is determined by 
dimensional considerations: 
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dK 
dy 



(3) 



where a is a dimensionless constant and (,{y) is the 'outer 
scale' whose physical meaning is the largest scale of tur- 
bulent fluctuations existing at distance y from the wall. 
We can define t{y) such that £{y) = y near the wall. Fur- 
ther from the wall l{y) saturates when coming close to 
the channel centerline. The full y dependence of (.{y) in 
three-dimensional channels was studied in with the 
final result l{y) ~ Lg {1 — exp [—A (1 + A/2)]} , where 
Ls = 0.3IIL and A = y(l - yl2L)/Ls. Note that the 
choice of Ls was made on the basis of 3D data, but we 
will keep the same value in 2D for lack of appropriate 
data. This will weakly affect the quantitative results but 
not at all the qualitative results below. 

In the vicinity of the wall, both in two and three di- 
mensions, we expect the flow to be differentiable, and 
moreover K{y) should start like This allows us 

to compute T>^ as —2vK{y)ly^. Since the energy pro- 
duction vanishes at the wall, this forces us to estimate 
the energy dissipation term near the wall as 



lim e{y) 



-2vK{y)ly^ 



(4) 



To model the energy dissipation away from the wall, 
we will ignore (for simplicity)the tensorial structure and 
write (up to a factor of unity) the kinetic energy dissipa- 
tion e ~ J/ J dkk^ IC{k) via the one dimensional turbulent 
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FIG. 1: Color online. Mean velocity profiles as a function of distance from the wall, in wall units, in three- and two-dimensional 
channels, for four values of the Reynolds number Re-r, shifted up by five units for clarity. The (grey) symbols in the left panel 
represent numerical simulations p^ . Note that the log-law (dashed lines) in 3D with an invariant von-Karman constant is 
increasing its range of validity whereas in 2D there is only an apparent log- law with a variable "constant" k (see insets). 



energy spectrum K,{k). Of course, this energy spectrum 
differs in two and three dimensions, resulting in a bifur- 
cation in the further development. 

3D case. In three dimensions the direct energy cas- 
cade results in the well known K41 spectrum /C3o(fc) — 
g 2/3 ^-5/3 [Y^^ rpj^g turbulent kinetic energy K{y) 



can be then estimated in the bulk as J^^^^dk JC3o{k) 
g2/3^2/3^y^ The upper limit was freely extended to infin- 
ity since the integral converges in the ultraviolet. Hence, 
e ~ 630 K'^^^/i. Joining up the estimates near and away 
from the wall we can write in three dimensions 

e{y) 0, 2 u [K{y)/l\y)] + 630 [K^/'iy) /^iv)] , (5) 

where near the wall £{y) y. 

2D case. In a two dimensional channel the inverse 
energy cascade does not play a role since the domi- 
nant driving is on a scale of the order of £{y). Away 
from the walls there exists a direct enstrophy cascade 
characterized by the Kraichnan [l^ energy spectrum 
K^u{k]y) ~ /32/3fc-3ln"^/3[H(y)], where (3 is the rate 
of enstrophy transfer. Therefore 

e{y) ^ V dkejC,^{k-y)^vl3'"H^^'^[£{y)/r^] 
Ji/t 

./3VMn^/3[^2(y)/3V3/,], (g) 



where rj ~ y/v/JP^ is the viscous (Kraichnan) length 
defined by the viscosity and enstrophy transfer rate. To 
eliminate /3 from these expressions we estimate K{ii) as 
JygdkJC2oik;y) ~ 13^/^1^ (y) for i{y) > ry. Using this in 
Eq. ^ we get 

e{y) = b2ui^Kr^l-a^^(v-H^/l<: + d^ . (7) 



Note that we have added a constant d in the argument 
of the logarithm which was only correct for large val- 
ues of i{y). With the regularizer we can go to the limit 
£(y) —^y—^0 without generating a spurious divergence. 
Requiring now that this equation agrees with the limit ^ 
we have to choose 62D = 2 and d ~ e. In other words. 



e{y) ^2,.[K{y)/i^{y)] (l{y)./K{y) / u ^ 
The resulting energy balance in 3D and 2D is 

WS+^(dlVK + A 



(8) 



dy \ J dy 

K . if3/2 



2v — \n^/^{v-^l^/K + c 



3D 
2D 



(9) 



Together with the momentum balance equation ([T]) we 
now have two equations relating our three objects S*, W 
and K. To complete the set of equations we employ a 
version of the Prandtl closure W{y) ~ h'T{y)S{y), which 
was carefully justified in (14i | for three-dimensional chan- 
nels, with the final result 



r„W^ceVKS, rw(2/) = [l + (^w/2/) 



6 1 1/6 



(10) 



where c and £^^1 are constants (£buf in wall units is 43 for 
3D current best fits). 

Having three equations for three unknown functions of 
y we can solve them given reasonable values of the param- 
eters. In three-dimensions we take the result of Ref. [l^. 
for the von-Karman constant K3D = (c^/fe)^/^. Experi- 
mentally in three-dimensions K3D ~ 0.415, and one is left 
with adjusting the constants a and b. We fix them using 
numerical simulation at the largest available Reynolds 
number, Re^ = 2003 [H, finding a « 0.218, b w 0.310. 
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FIG. 2: Color online. Left and middle panels: Profiles of the turbulent kinetic energy K'^ — K/{p'L) in 3D and 2D channels 
respectively. Symbols in the left panel are from numerical simulations Right panel: Reynolds shear stress = W/{p'L) 
as a function of y/L in three (solid lines) and two-dimensions (dashed lines), respectively, for four different values of Rct. The 
symbols are numerical simulations 13,]. The inset in the left panel compares the profiles of K in three- and two-dimensions at 
Re,- = 2003 normalized to their maximal values. 



Not having independent data in two-dimensions we take 
there the same values for a, and Lg. It was suggested 
in Ref. [uj that k^d « 0.2 for Re^ « 10^. While we do 
not agree with this reference on the existence of a power 
law in 2D, there is a small range of y where such a law 
can be fitted also in our results, and we use this number 
to adjust the value of c to c w 0.047. This completes the 
choice of parameters in two an three dimensions. 

The theoretical predictions for the mean velocity pro- 
files are shown in Fig. [1] where we have used the wall- 
coordinates RBt = Ly/p'L/iy, y+ = j/Re^/L, = 
Y I yJp'L. Note the good agreement between predic- 
tions and data in three-dimensions throughout the en- 
tire channel, including the viscous, buffer, log-law and 
wake regions. This underlines the quality of the model 
employed here. We note that in 3D the log-law region 
(shown as dashed line) increases with Rcr. Indeed, the 
present model is very similar to the one proposed in 14 1 
in which the von-Karman log-law k. « 0.415 was ob- 
tained. The model also captures the "wake" feature of 
the velocity profiles. In contrast, the mean velocity pro- 
files for two-dimensional channel in the right panel reveal 
a very limited log-law region (if any) as well as an increase 
in the apparent von-Karman constant upon increasing 
Rer, see inset in Fig. [T] right panel. The "wake" feature 
is also absent, instead the curves bend down. 

Even a larger qualitative difference is seen in the tur- 
bulent kinetic energy profiles in wall units, as seen in 
Fig. [21 left and middle panels. We first note the order 
of magnitude difference in the value of the kinetic en- 
ergy in favor of the two-dimensional channel. Second, 
the 3D profiles are almost Rei--independent, saturating 
at Rcr — *■ c». In contradistinction the 2D profiles in- 
crease almost linearly with In(ReT-). The third difference 
is the decline of the 3D with i/+ outside the buffer 
layer, reaching a pronounced minimum at the centerline. 
In the two-dimensional case becomes almost in- 



dependent. Interestingly enough, the highly different be- 
havior of is not mirrored in the profiles of which 
are shown in the right panel of Fig. [51 

In order to rationalize all the differences presented 
above, we will discuss the energy balances between pro- 
duction, diffusion and dissipation which are shown in the 
left panel of Fig. [3l In order not to get mixed up be- 
tween diffusion and dissipation we will consider the total 
(integral over the channel) energy balance from which the 
diffusion term disappears exactly. The integral of energy 
production Ep = WSdy is dominated by the bulk re- 
gion y > y* = y^p' L/v where y^ ~ 20. In this region we 
can take W ^ p'L and S{y) ~ y/p'L/{Ki'y) where k is 
the von-Karman constant in three-dimensions and a very 
weakly Re^ dependent number in two-dimensions. Inte- 
grating between and L we end up with the estimate 
[(p'L)3/2/«] ln(L/y,) or 



[(p'L)3/V«] ln(Re./y^) 



(11) 



On the other hand the total dissipation Ejj differs in 
three- and two-dimensions. In 3D we estimate the 
main term from the bulk of the channel as Ed « 
b K^/'^{y)/y w K^/'^\n(Rer/yt), where K is some 
typical value of K(y) in the bulk. Comparing with 
Eq. pT|) we see that indeed K ~p'L and therefore Re,-- 
independent, cf. Fig. [21 left panel. This is not the case in 
two-dimensions. According to Eq. ([9]) the dissipation in 
2D is roughly constant up to y ~ y ~ 10 p'L/v, whereas 
for y > y it decreases roughly like y^^. Moreover, the 
integral over the dissipation is roughly equal in these two 
regions. We can therefore estimate -Ed ~ -ii^K J~ dy/y^ 
where the weak logarithm is neglected. This integral 
is dominated by its lower limit and therefore in two- 
dimensions -Ed ~ AvK/y ^ K^Jp'L/y^ . Comparing 
with Eq. (fTT|) we sec that 



{y+/K) In (Re^/y+) - A In (Re^) + B , (12) 
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FIG. 3: Color online. Wall units. Left panel: The energy balance between production, diffusion and (minus) dissipation for 3D 
(solid lines) and 2D (broken lines) channels at Re^ = 2003. Middle panel: The Rer dependence of typical values of the kinetic 
energy in 2D channels. Note the agreement with the estimate in Eq. (|12p . Right panel: Averaged turbulent kinetic energy over 



the channel half- width L: K^, 



L ^ K{y) dy. Inset: the quick saturation of A'^gj^go in 3D channels. 



with A and B being constants. This explains nicely the 
results of Fig. [2l middle panel, where we see that mul- 
tiplying the values of Re,- by a factor of five results in 
equidistant curves of in linear scale, and this is why 
K can be taken as the asymptotic (almost constant) value 
of K. To support these estimates by numerics we show 
in Fig. [3l middle panel, the actual functional dependence 
of i^max ^^-d at various position in the channel on 
Rcr. Equation is very well supported by these data. 

Finally, we want to explain the profiles in Fig. [TJ The 
log-law in three-dimensions can easily be derived by tak- 
ing the dominant terms in Eqs. ([T]), ([9|) and (fTO|) (in 
the limits Re,- — > oo and p' L / v <^ y <^ L)\ W — p'L, 
WS = hK^I'^/y and W = cy^KS. Solving these for S 
we find S{y) oc l/y i.e log profile. This simple solution 
is lost in 2D because of the absence of the direct cascade 
that lead to the term bK^^'^ /y. Nevertheless we see that 
the role of the direct cascade in two-dimensions is mim- 
icked by the (negative of the) diffusive term which is also 
Re^. independent and leading to a similar estimate if we 
replace dK/dy by K/y. For that reason we still have a 
remnant of a log-law also in 2D, but with a slope that 
changes with Rcr due to the dependence of K on Re^. 

In conclusion, we show that the Reynolds stress profiles 
in 2D and 3D channels look similar. The velocity profiles 
in 3D are truly represented by a log-law, but in 2D they 
are only apparently varying according to a log-law, the 
"constant" is changing as a function of Re^, even though 
not vary rapidly. But the major change between two and 
three dimensions is in the kinetic energy profile. The 
two-dimensional channel is much more energetic, with 
the mean kinetic energy increasing like ln(Re7-). This 
is the main result of the loss of the direct energy cas- 
cade. Nevertheless even without this cascade the energy 
dissipation exceeds its laminar value and it changes the 
dependence of the kinetic energy on Re,- from 0(Rer) 
as in laminar flows, to ©(InRcr) which is one of the in- 
teresting predictions offered above. To make this point 
clear we show in the right panel of Fig. 3 the difference 



between the averaged kinetic energy in the channel in 
2D and 3D as a function of Rcr. The former is roughly 
linear in In(Rer) and the the latter saturates quickly to 
a constant of about 2. It should be tempting to test in 
experiments some of these predictions, even if the ideal 
2D channel model is not easy to achieve. 
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